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О стойкости методов встраивания 
цифровой информации к атакам 
сглаживающими фильтрами 


Рассмотрены общие свойства линейных сглаживающих фильтров в области Фурье и в пространственной 
области. Основное внимание уделено сохранению параметров квазилинейных и квазиплоских сигналов. 
Доказаны теоремы о сохранении близости квазилинейного сигнала к линейной функции и квазиплоского 
сигнала к функции плоскости. Получены верхние оценки разности оценок коэффициентов наклона 
линейной функции к оси абсцисс до и после работы сглаживающего фильтра. 


Введение 


В цифровой обработке изображений (ЦОС) класс операций, который соответствую- 
щим образом объединяет пиксели малой окрестности исходного изображения и форми- 
рует новое изображение, называется классом операций над соседними элементами [1]. 
Эти операции применяются при низкоуровневой обработке изображений и называются 
фильтрами. Дискретный оператор Н формирования окрестностей изображения ото- 
бражает матрицу размерностью МХМ в саму себя операцией 

С, =Н(С )5 Ут ‚п ] ЕЁ, (1) 


т-т,п’-п 
где Г — дискретное множество точек, называемое окном или маской; бы О т 
матрица пикселей изображения размерностью МхМ до и после применения фильтра. 
Элементарная комбинация пикселей в окне задается умножением величины 
каждого пикселя в пределах маски фильтра на соответствующий весовой множитель 
маски и сложением полученных произведений по всей маске. Результат присваивается 


значению центрального пикселя: 


р с) 


т=фки=-г 


где М =М= 2+ 1, и — весовой коэффициент для элемента изображения ый 


попадающего в маску с центром в точке с координатами (т,п). 
Операция (2) эквивалентна операции дискретной свертки 
М> м-— 
Эти = зе Е | 
т =0п =0 
Развитие цифровой обработки сигналов основывалось на использовании дискрет- 
ного преобразования Фурье периодических функций. При применении преобразования 
Фурье к конечному дискретному упорядоченному множеству, представляющему 
сигнал, это множество отождествляется с периодом бесконечной периодической функ- 
ции. В одномерном случае бесконечное число раз справа и слева к исходной конеч- 
ной последовательности добавляется та же самая последовательность. В результате полу- 


«Штучний 1нтелект» 372009 587 


Никитенко Л.Л. 


тон 


чается бесконечная периодическая последовательность, один период которой вклю- 
чает всю исходную конечную последовательность. 

Дискретное преобразование Фурье отображает матрицу изображения размерностью 
Мхм в комплекснозначную матрицу той же размерности: 


| >> 28 и, ий 
Зиму ММ 8т‚п @Хр М р де” 


т’=0 п’ =0 


где &„ — элемент матрицы изображения в области Фурье. В соответствии с теоремой о 


свертке в области Фурье свертка сводится к умножению: 
С*Н <> ММСН. 

Коэффициент ММ компенсирует усреднение по матрице изображения. Преобразо- 
вание Фурье от маски свертки Н называется передаточной функцией фильтра. В общем 
случае применение фильтров меняет и амплитуду, и фазу элементов изображения. 

Целью данной работы является поиск сохраняющихся или мало изменяющихся 


параметров цифровых сигналов при обработке линейными сглаживающими фильтрами 
для создания стеганосистем, стойких к атакам сглаживания. 


Общие свойства сглаживающих фильтров 


При низкоуровневой обработке изображений с целью уменьшения влияния шума, 
неравномерной освещенности, неоднородного фона и т.д. было создано множество раз- 
нообразных фильтров, которые называются усредняющими или сглаживающими филь- 
трами [1]. В области Фурье передаточная функция линейного сглаживающего фильтра 
имеет четыре свойства, обеспечивающие как можно лучшее сохранение качества исход- 
ного изображения при уменьшении влияния высокочастотных элементов: 

1) действительная, не вносятся изменения в фазу элементов изображения; 


2) сохраняющая среднее значение, т.е. Про =1; 


3) изотропная, сглаживание одинаково во всех направлениях; 

4) убывающая с ростом волнового числа (уменьшает влияние высокочастотных 
элементов). 

Маска размером (2^ +1)» (2 +1) с четной горизонтальной и вертикальной 


симметрией приводит к передаточной функции: 


у К 
#(К) = йо + в соз(и'лЁ! ) + 2 йо соз(т ло) + 
п’ п’ = 
ЗИ , , (3) 
+4 >. ьз п’, с0$(и лА1 ) с0$(т ло ), 


т =Ш =1 
где К, и К, — соответствующие волновые числа по горизонтали и вертикали. 


Если сглаживающий фильтр применяется после того, как изображение было по- 
делено на непересекающиеся блоки с размерами, совпадающими с размерами маски 
фильтра, то фазы элементов в области Фурье и нулевой элемент (среднее значение) 
останутся неизменными. Значит, в этом случае дополнительная информация, встроен- 
ная в нулевой элемент и в фазы элементов блока, не будет подвергаться искажению. 

Можно сделать вывод, что методы, встраивающие дополнительную информацию 
в нулевой элемент и в фазы элементов блоков изображения в области Фурье, будут 
стойкими к атакам линейными сглаживающими фильтрами, если 
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1) размеры маски фильтра совпадают с размерами блоков, на которые разделялось 
изображение при встраивании дополнительной информации; 

2) в пространственной области разделение изображения на блоки при применении 
сглаживающих фильтров совпадает с тем разделением на блоки, которое выполнялось 
при встраивании дополнительной информации. 

Понятно, что перечисленные условия, обеспечивающие стойкость методов встраи- 
вания к подобным атакам, на практике трудно осуществимы. Поэтому целесообразно 
вернуться в пространственную область и рассмотреть работу сглаживающих фильтров с 
самого начала. 

В пространственной области для сглаживающих фильтров выполняется условие 


о =1, из которого и получается условие Жо=1 в области Фурье. Это 


т п 


условие, как правило, выполняется и для нелинейных сглаживающих фильтров, 
например, для медианного фильтра [1]. Кроме того, остается обязательным и условие 
изотропности фильтра. В реальных линейных фильтрах условие изотропности всегда 
выполняется для вертикального и горизонтального направлений, а по другим направле- 
ниям к этому условию стараются максимально приблизиться. Далее рассматриваются 
только линейные сглаживающие фильтры. 

Четные сглаживающие фильтры [1], в отличие от нечетных, не сдвигают элементы 
исходного изображения с их позиций в цифровой матрице сигнала. В данной работе 
будут рассматриваться только четные фильтры, применяемые для одно- и двумерных 
сигналов. В общем случае влияние линейного четного сглаживающего фильтра для 
одномерного сигнала имеет вид: 


8. а п 8, ы >.*.. 2, ы а ) > (4) 


п’ = 


а для двумерного сигнала 


Я =. По 8 ии ы 2 Вия + ты ) ы йа и + РН т 
п’ =1 т =1 
(5) 


Г Г 
т уз 2% п вы ы и а ы т ). 


т =1и = 

Пусть во время работы сглаживающего фильтра происходит последовательная 
обработка всех пикселей исходного сигнала. Если сигнал описывается линейной 
функцией (в одномерном случае) или функцией плоскости (в двумерном случае), 
тогда сглаживающий линейный фильтр не изменит его значений. Кроме того, 
остаются неизменными и параметры (коэффициенты наклона к осям абсцисс и 
свободный член) линейных (плоских для двумерного сигнала) функций, описывающих 
значения величин пикселей в блоке. 


Сглаживание квазилинейных сигналов 


Пусть последовательность величин пикселей в одномерном сигнале 2=(®,&,..2\, |) 
не является точной дискретной линейной функцией / = АТ+С, но близка к ней. 


Формально близость последовательности к линейной функции будем понимать в том 
смысле, что максимальное отклонение элементов последовательности отличается от 
соответствующего значения функции не более чем на заданную величину Д: 


|=, (6) 
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Такой сигнал будем называть квазилинейным. Применим линейный сглаживающий 
фильтр к последовательности элементов одномерного квазилинейного дискретного 
сигнала. Очевидно, что последовательность на выходе сглаживающего фильтра тоже 
будет квазилинейной. Тем не менее, этот факт требует четкого математического 
обоснования. Кроме того, необходимо оценить, насколько близкой к линейной 
функции будет выходная последовательность по сравнению с близостью к этой 
функции входной последовательности. 

Теорема 1. Сглаживающий фильтр сохраняет близость квазилинейного сигнала 
к линейной функции, т.е. 


1 
|=, - 15 А> |8!-Л|<А, (7) 
где а, ©; — элементы последовательности сигнала на входе и выходе сглаживающего 
фильтра соответственно. 

Доказательство. Дополним последовательность © = ($1, 22....2м) справа и слева 
значениями функции /= К; + С. Введем обозначение Д; = ©; и применим произвольный 
сглаживающий фильтр к полученной последовательности. 

и 
= Г/) + У» ' "+ ") = 
&п = т Вы) 


- п-+п 
п =1 


Ро ре тк г а. м — 
п’ 
= (Кп+С+А,) +». й,(К(п-п)+К(и+п')+2С+А_,+А,,)= 
п’ = 
= (Ки +С+)+>`й,(2Кп+2С) + А+ й,(А,„+А,,)= 


п =1 п =1 


= +И А, +У'^, (А_.+А,)=Х, +А, . 
"= 


й 
Оценим величину Д, ‚ используя свойство сглаживающих фильтров Ув =| и 
т=-г 


неравенство (6): 


А, =ДА, +УЙ. (А о +А .)<АА+2У Й.А=АВ +2 й.)>А. 
п 0 п п п-п п+п 0 - п 0 - п 


п’ = п =1 п =1 


Теорема доказана. 

В формулировке и доказательстве теоремы 1 использованы известные параметры 
(коэффициент наклона К и свободный член С’) линейной функции, к которой близка 
последовательность © = (2,,8.,...2».). На практике эти параметры чаще всего не 


известны. Более того, получателю часто не известна сама последовательность 
8 = (2,8,...8 у), ему известна последовательность на выходе сглаживающего фильтра 


2 = (2%,91,...8 м1). Значит, следующим шагом в исследовании влияния сглаживаю- 
щих фильтров на квазилинейную последовательность должна быть оценка коэффициента 
наклона линейной функции (далее коэффициент наклона) до и после работы сглаживаю- 
щего фильтра и их сравнение. 

Оценим свободный член линейной функции по формуле 


= в М-1 М-—1 
С = К — К й 8 
м2. 8, Е - (8) 
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а коэффициент наклона с помощью метода наименьших квадратов 
УК б- 4.) ва. (9) 
п=0 
С учетом (8) выражение (9) можно переписать в виде 
У (Киа. К = 


Продифференцируем левую часть (10) по К и приравняем производную нулю, 
чтобы получить оценку К в явном виде: 


ны Иа М1 
2 (К(и- и а, 
п=0 


2,) > ша. (10) 


КУ 7: =>. в -^ 5. (т 
Поскольку [2] 
с — т м у й щи 


(12) 


= 
[ 
- 


то 


№-—1 = №-—1 =] —] 
Е ос в 
п=0 м. п=0 2 2 


_ мм-@м-1 и А-В, 0 _ 
6 2 4 

_ Мм- 
12 

_ ММ-П(У+ _ ММ? -1 

* 12 Е И. 


Правую часть уравнения (11) тоже можно упростить 


2.4в.- во = Ув. -^—- Уи -^—= 


(22№М-1)-6(№М-1)+3(№-1) = 
(13) 


—1 М-—1 
_1 =] 1 
№ [9 и т „мм 


(14) 


Подставляя (13) и (14) в (11), получаем окончательное выражение для оценки 
коэффициента наклона линейной функции, к которой близка последовательность 


& = (80, 81›--.8 м): 
- 12 т М(М-1) 
К=— пе. 15 
(в =) (15) 
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Выражение (15) можно представить в другом виде, который более удобен для 
последующих исследований: 


г _ \; ее р 
= в. | ыы 


Оценка К вычисляется по (15) и (16), но вместо ©, подставляется ©’. 


Найдем погрешность А, полученной оценки К’ по отношению к оценке К : 
М- 
(М-| 
ды в] 
№ т = 2 


Под знаком суммы с правой стороны (17) в качестве сомножителя стоит 
величина (2, — 2, ), которую стоит рассмотреть отдельно: 


(17) 


(2, -2,) = (Кп+С+А, -(Кп+С+А,)) ЕД, -А,. 


Подставим в (17) полученное значение (2, -©,„)и сделаем предположение, что 


М - четное число. (В цифровой стеганографии практически всегда выбираются 
последовательности с длиной 27, где р — целое число.) 


г 4-9] 


К МОП 


св (-°] 


Е 2 
М. 155 (М-1) 10. || _ (м-9)_ 
< т 2.4 + ы - "| ат 24 + 5}, —- = 


244 | ((М-0 А [М _^)_(мМ- _ 
а, 2 "В+ } 2 
___ 24А и 0+ 24Л [АНА 1М№\ _ 
ММ? -12 2 2 2 Ма Пе р 


= ты [м-0-(5-9+[(%+)-5- р 
А4М(М№? —1) 2 2 (М -1 


Полученный результат сформулируем в виде леммы. 
Лемма 1. Если параметры линейной функции /, = КТ + С ‚ описывающей дискрет- 


1 


ную квазилинейную последовательность © = (85, 81»... и_1), оцениваются по формулам: 


Ув. и 1) 
= = | 
а м 
С= К 

м2. 2 


НИИ 
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и в, - (Кп+С) <А, тогда 


а (18) 
(М -1 
Из (18) видно, что при М№М->ю АД, ->0. В частности, при М№М=8 
ль 0,764. ‚ апри №М=16 ДА ды де 0,38 д. 
63 255 


Сглаживание квазиплоских сигналов 


По аналогии с квазилинейными сигналами введем понятие квазиплоских дис- 
кретных сигналов, т.е. двумерных дискретных сигналов, в которых значения пикселей 
2„„ Удалены не более чем на величину Д от соответствующих значений плоской 


дискретной функции /„,: 
/„ =Кт+К,п+С, (19) 


[т — /ш| ЗА. (20) 


Теорема 2. Сглаживающий фильтр сохраняет близость квазиплоского сигнала к 
плоской дискретной функции, т.е. 


<А, (21) 


Я 


|2 т — У, г тш <А >> [= _— 


где 2„,©, - элементы сигнала на входе и выходе сглаживающего фильтра 


соответственно. 
Доказательство. Дополним прилежащие области к множеству элементов сигнала 
С={е „} значениями функции /„, =К\т+К,п+С. Введем обозначение 


А = 8 тп В И (22) 


и применим произвольный сглаживающий фильтр Н к полученному множеству 
элементов сигнала: 


Г № Ни Е А } = и Е А 2 (23) 


где в соответствии с (5) 


„м = Яо и +, ‚(А т‚п-п’ А ди ‚+ >, о(Д т-т ‚п АЕ 


п’ = т =1 


в Ш 


+ 2 уз п, п’ (А ы В т ли + В ). 


Верхнюю оценку отклонения элементов сглаженного сигнала от соответствующих 
значений плоской дискретной функции можно получить, используя неравенство 


| 


< А, которое тождественно неравенству (20): 


Ат < Йо А +, 2А+Ул, АУУ, ‚4А = лу Ул, 


т =! п'=Ёи = п=фит=-ф 


тп 


Теорема доказана. 
Теорема 1 является частным случаем теоремы 2, но она может иметь большое зна- 
чение в теории использования аудиосигналов, поэтому ее доказательство было приве- 
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дено в полном объеме. Кроме того, последовательность действий в доказательстве тео- 
ремы 1 позволяет сократить последовательность действий при доказательстве теоремы 2 
и облегчает понимание смысла последовательности тождеств (23). 
Оценки коэффициентов наклона А, и К, плоской функции и свободного члена 
С получим с помощью метода наименьших квадратов: 
М-1 №М-1 
У» (Кт+Ки+С-8„) > тм. (24) 
т=0 п=0 
Чтобы оценить С, продифференцируем (24) по С и приравняем производную нулю: 
МЬ-—1 М-1 


УУ(Кт+Ки+С-8„)=0, 
т=0 п=0 
— 1 МГ №М-—1 М-— 1 
бету аи ++). 65) 
т=0 п=0 


Иначе это равенство можно записать, используя среднее значение элементов 
двумерного сигнала &„: 


Сакко. 


(26) 


Подставим оценку С в (24) и снова продифференцируем (24) по К, и К.. 
Приравняем производные нулю и решим полученную систему равенств, чтобы найти 
оценки А; и К, 


У (к за и" Е в, мк ия ==} а 


Ни ` . (27) 
т=0и=0 о 2 - Е 
Там, где это возможно, вынесем за знак суммы постоянные величины: 
М-1 М-ЕМ-1 к 
КМУ т- ео ] +К У (»- т ии °)- 
т=0 т=0п=0 2 
М-1м- М1 
сер 
т=0и=0 
М- М-1М-1 ки 
км» [„- ей ] + +К РУ [#9 ии 2) 
т=0 т=0п=0 р 
М-ЕМ- 
(и-1 | 
= — и- 
т и 2 
Чтобы остить полученную систему, воспользуемся равенствами (12) и 
М-1 
(М-П : 
авенством ЕТО 
: > 2 
№2 м2 —1) ЗЫ М-1 
К 12 0 у [ив 
т=0п=0 : 
М2[№2 1] _ У М-1 
вы уз Хо ("- 2 
т=0п=0 
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О стойкости методов встраивания цифровой информации... 


тон 


ИЛИ 


(28) 


Из представления -. видно, что оценки коэффициентов наклона плоской 
функции напоминают оценку коэффициента наклона в линейном случае, но 
добавляется усреднение по второй координате. 


В общем виде, как и следовало ожидать, выражения для К, и К, совпадают, 
поэтому достаточно сравнить оценки одного коэффициента, например К,, до и 


после работы сглаживающего фильтра. Оценивание ЛД, проведем по аналогии с 
исследованиями в случае одномерного сигнала: 


_ м ММ 
Акт = К |= 8 =) [= и 
лв 
М-Т |2 №—1 
12 у —1 у М-1 
=— _ + ( = Дт-—— || 
2.2 > > бы 5 Е > тп 5 
№ (М -1 п=0 ||т=0 И 2 
2 
Поскольку 


8 тн т рт = т А их НХ 5 А) = у 7% А < 2А ? 
то 
ИЗ 


Е И > 


(№ =П п=0 т=0 


мА 
ео 


ея и 
Выражение для верхней оценки Л к Для двумерного сигнала совпадает с 
выражением для верхней оценки погрешность АД; для одномерного сигнала. 
Результат сформулируем в виде леммы. 
Лемма 2. Если параметры плоской дискретной функции / 


тп 


=Кт+К,п+С, опи- 
сывающей дискретное квазиплоское множество С’, = {2 „ !, оцениваются по формулам 


МЬ-1М-—1 
ы м р 
К ЗаяОи= =) 
=» 
_ 12 МЬ-М-1 ть В ? 
КЕ Е. _ 
2 г (№? ы 1) 2. 2.втя (п ) 


а що, 


т=0п=0 
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и |» -(Кт+К.п+ С) <А, тогда различие Ан, #= {1:2} между оценками коэф- 


фициента К, до и после обработки сглаживающим фильтром ограниченно согласно вы- 
ражению: 


Е (29) 


Выводы 


В работе изучалась ЦОС сглаживающими фильтрами элементов и параметров 
сигналов в случаях одно- и двумерных сигналов. Были рассмотрены общие свойства 
линейных сглаживающих фильтров в области преобразования Фурье. Основное вни- 
мание уделялось фильтрам с четной горизонтальной и вертикальной симметрией, по- 
скольку их применение не вызывает смещения сигналов в пространственной области. 

В пространственной области для сигналов, которые можно считать линейными в 
случае одномерных сигналов или плоскими в случае двумерных, были доказаны теоремы 
о сохранении параметров линейности (плоскости) при применении линейных сглаживаю- 
щих фильтров. Кроме того, в обоих случаях были получены верхние оценки отклонения 
оценок коэффициентов наклона, рассчитанных после применения сглаживающего фильтра, 
от оценок, рассчитанных до применения сглаживающего фильтра. Оценивание коэффи- 
циентов наклона проводилось методом наименьших квадратов. 

Исследования позволяют сделать вывод, что если сглаживание будет произво- 
диться в области Фурье с заранее известными блоками, то для встраивания секретной 
информации можно использовать нулевые элементы блока и фазы остальных элементов. 
При сглаживании в пространственной области возможно встраивание дополни- 
тельной информации в коэффициенты наклона линейных функций, описывающих 
квазилинейный сигнал, и в коэффициенты наклона функций плоскости, описыва- 
ющих квазиплоский сигнал. При этом следует учитывать погрешности, возникающие 
при определении коэффициентов наклона получателем после применения сглажи- 
вающих фильтров. 
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Про стиксть методв вкраплення цифрово! 1нформацИ щодо атак згладжуючими ф!льтрами 
Розглядалися загальн! властивост! лийних згладжуючих фильтрв в област! Фур’е та у просторовйй област. 
Основна увага придлена збереженню параметрав квазилийних 1 квазшлоских сигнал. Доведено теореми 
про збереження близькост! квазиийного сигналу до лиийно! функцй та квазшлоского сигналу до функцИ 
площини. Отримано верхн! ощнки рзнищ! ощнок коефищенпв нахилу до та шсля роботи згладжуючого фильтра. 


[Г.. МийепКо 

Оп Фе Эа у о Фе Нше Ме@од$ о? е Попа Мез5асе уу Вебрес 0 ЭЗтоо то ЕШег$ 

ТБе рарег деуое4 ю Фе Чеиа| $1юпа| ргосеззте плефо4$ гезеагсЬ Бу зтооте Н\егз. Сепега| ргорегаез оЁ пе 
зтоофте ННегз ууеге ехаттед ш Фе Еоипег доташ ап ш Фе зрана| 4оташт. СепегаПу, сБагасет$с$ 
тезегуаноп оЁ фе диазШпеаг 91па[$ ап4 аиазр]апаг па! уме аНепаед. ТВе Феогепа$ абойй гезегуаноп оЁ фе 
ЧлазШпеаг $1опа] ргохнайу 0 фе Ппеаг КлпсНоп ап4 доазр]апаг 51а] ргохнайу ю фе рапе вшсйоп \еге ргоуе4. 
Тре иррег Боип4$ оЁ Фе Ч!егепсез Бебмееп Ппеаг КпсНоп шсПпайоп соеЁйстеп( езбтланопз ю фе абзс1зза ахез 
БеГоге ап4 айег фе злооф те В ег асноп \уеге гезиНед. 
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